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The article continues the author's publication in [Mech. Tverd. Tela, No. 34, 2004], in which the generalizations of 
the Appelrot classes of the Kowalevski top motions are found for the case of the double force field. We consider the analogue 
of the 4th Appelrot class. The trajectories of this family fill the surface which is four-dimensional in the neighborhood of 
its generic points. The complete system of two integrals is pointed out. For these integrals the bifurcation diagram is 
established and the admissible region for the corresponding constants is found. 



1 Introduction 

Consider a rigid body with a fixed point O. Let the principal moments of inertia at O satisfy the ratio 
2:2:1. Suppose that the moment of external forces with respect to O has the form 

ei X a + e2 X /3, 

where the vectors ei, 62 are fixed in the body and parallel to the equatorial plane of the inertia ellipsoid, 
and a,/3 are the vectors constant in the inertial space. 

It is shown in the work [5] that without loss of generality one can consider ei,e2 to form the 
orthonormal pair (in particular, to be the principal inertia unit vectors), and a, (3 to be mutually orthog- 
onal. Let 63 = ei X 62. Choose 0616263 as the moving frame. Denote by U3 the angular velocity vector. 
In the dimension less variables the rotation of the body is described by the Euler-Poisson equations 

2loi — aj2^3 + Pz, 2a;2 — —^1^3 — 013, ^3 = Q!2 — 

(1) 

ai = a2UJ3 - Pi = ^2t^3 - /^3W2 (123) . 

The (123) symbol means that the remaining equations of the Poisson group are obtained by the cyclic 
substitution of the indexes. 

The phase space of system ([T]) is defined in R^(u;, a,f3) by the geometric integrals 

al + al + al = a\ f3f + (3i + f3i = b\ ai(3i + a2p2 + a^Ps = 0- (2) 

We suppose that 

a > 6 > 0. (3) 

Then system ([1]), ([2]) does not have any cyclic integrals and is not reducible, by the standard procedure, 
to a Hamiltonian system with two degrees of freedom. Nevertheless, it is completely integrable due to 
the existence of the first integrals in involution 

H ^ Ljj + uj'^ + ^ujj - {ai + (32), 

K ^ {ujl^Lul + ai- P2f + (2cjiW2 + a2+ /3l)^ 

(4) 

G = {aiuji + 0:2^2 + ^ctsujs)'^ + {Piuji + ^2(^2 + 5/^3^3)^+ 
-I- ^3(71^1 + 72W2 + 573^3) - "1^^ - P2a^ 
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(here 7^ stand for the components of the vector cy. x (3 immovable in space). 

The integral K was first shown by O.I. Bogoyavlensky [2], and the integral G (in more general form 
for a gyrostat) was found by A.G. Reyman and M.A. Semenov-Tian-Shansky [4]. 

In the work [S] the set of critical points is found for the integral map 

H X K X G : P^' ^R.^. (5) 

It is shown that this set is the union of three sets 01, D, which are almost everywhere the smooth 
four-dimensional submanifolds in and in the neighborhood of the generic points are defined by two 
invariant relations. 

The first critical set DJfl was found in the work [5]. It coincides with the zero level of the integral 
K and generalizes the l*** Appelrot class in the Kowalevski problem [T]. The phase topology of the 
dynamical system induced on 971 was studied by D.B. Zotev 0. 

The motions on the manifold found in [6] are investigated in [7l [8] . It is shown that this family of 
motions is the generalization of the so-called especially remarkable motions of the 2"*^ and 3'''^ Appelrot 
classes. The bifurcation diagram for the pair of almost everywhere independent first integrals on *Jt is 
constructed in the equations on ^ are separated, the bifurcations of the Liouville tori are studied in 

The present work is devoted to the investigation of some properties of system ([T]) restricted to the 
invariant subset D given in P^ by the system of the invariant relations [5] 



where 



Ri - 0, i?2 = 0, (6) 

Ri = (a3W2 - P3^i)^3 - 2/3iCjJ + 2{ai - (32)ujiuj2 + 2a2a^|, 
i?2 = {a^^i + I33UJ2)ujI + [al + /3| + 2aiujl + 2(a2 + /3i)t^iW2 + 2/32t^i]t^3+ 
-I- 2a3[(ai - /32)wi + (^2 + /3i)w2] + 2/33[(a2 -|- /3i)wi - (ai - I32)lo2]- 

2 Partial integrals 

Note that at the points 

UJi= tjJ2 = 0, Q!3 = /33 = (7) 

equations (O are dependent on D. If we assume that equations ([7]) hold during some time interval (and 
then identically in t along the whole trajectory), then we come to the family of pendulum type motions 

a = a(ei cos — 62 sin ^^), /3 ±fe(ei sin6' + 62 cos6'), a x /3 = ±0663, 

(8) 

u} = ee3, 6* = -(a± 6)sin6' 
noticed in the work [5^. The constants of integrals (H]) at such trajectories satisfy one of the following: 

g = abh, k^{a-bf, h^-{a + b) (9) 

or 

g = -abh, k^{a + bf, h^~{a-b). (10) 

Denote by the set of points belonging to trajectories Let D* = D\^. 

Recall that in the classical case of S. Kowalevski {f3 — 0) there exists the area integral. Traditionally 
it is represented with the one-half multiplier 

L = -Iu}-a. (11) 
2 ^ ' 
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Then provided that /3 — the integral G turns into L^. 

Let £ be the constant of integral According to G.G. Appelrot's classification the 4**^ class of 

the especially remarkable motions is defined by the following conditions: 

1) the second polynomial of Kowalevski has a multiple root, one of the Kowalevski variables 
remains constant and equal to the multiple root s of the corresponding Euler resolvent defined as 
ip{s) = s(s - h)^ + (a2 - k)s - 2^2. 

ip{s) = 0, ip'is) = 0; (12) 

2) the first two components of the angular velocity are constant and equal to 

e , , 

wi = --, u)i = 0. (13) 
s 

The next statement establishes the analogue of conditions (fT3|) for the generalized top. 

Theorem 1. For any trajectory in the set O* the values and are equal to each other and 

Iw-ei \ijj-e2 

constant. 
Proof. Denote 

M = lu;, 

Ml = Iu)-ei = 2wi, M2 = Iu)-e2 = 2uj2, 
Ma = lij-a = 2q:icji + 2a2UJ2 + aa^s, — li^-fH = 2/3ia;i + 2/32^^2 + /Sa^s- 

The first equation ^ becomes 

^-^.0. (14) 
Ml M2 

Introduce the function 

MaMi + MpM2 



^ Ml + Af| 



Its derivative in virtue of equations ([T]) is 
dS 1 



dt 4(M2 + Af I 



[{M( + M|)w3 + 4a3Mi + AI3iM2]{MpMi - MaM2). 



Then (|14p implies that the right hand part vanishes identically. Therefore, S is the partial integral on 
the set D* . Let s stand for the constant of this integral, 

M^Mi + Mf,M2 _ 

Ml + M| ^ > 

From we obtain that A/q = — sAfi, Af^ = —SM2 with the constant value s. □ 

Remark 1. In virtue of condition (fT4l) the function S* can be also written in the form 



l,Ma Mq^ , , 

2(a^+a|)- 

Remark 2. Note the interesting geometric feature of the kinetic momentum vector motion on the 
trajectories considered. Introduce the immovable orthonormal basis in the Oct^-plane 

a j3 

1^1 = — , V2 = -. 

a 
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Let mi = M'l^i, m2 = M'i/2. Then Ma = ami, Mp = bm2, and condition ((T4|) yields 

Af2 m2 
Ml a mi ' 

or 

tgi?= - tg7?o, 
a 

where i?o are the polar angles of the projections of the vector M, respectively, onto the equatorial plane 
of the body and onto the plane of the direction vectors of the forces fields. 

Theorem 2. On the set D system ([7]) has the partial integral 

T = {a^oji + l33Uj2)oJ3 + 2aiujf + 2{a2 + /3i)uJiti'2 + 2/32ti'|— 
-2(ai/32-a2/3i) + a2 + fo2. 

Proof. The derivative of (fT7)) in virtue of system ([1]) is 



(17) 



dT 1 

and vanishes identically on D. □ 
Denote by r the constant of the integral T. 

In the work j5] equations ([6]) are obtained from the condition that the function with Lagrange's 
multipliers s, t 

2G+{t-p^)H + sK 

has a critical point. Comparing (|16p . (jl7p with the expressions for s, t in |5], we see that these multipliers 
are the constants of the above given integrals S, T. 

According to ^ introduce the positive parameters p, r as follows 

2 2 I 7,2 2 2 J 2 

p=a+o, r = a — . 

Let h,k,g be the constants of the general integrals Then equations (50) of the work [S] give the 
following relations on the set D*, 

These relations also can be considered as the parametric equations of the sheet of the bifurcation diagram 
of the map ([5]). Eliminating of r leads to the equations 

V'(s) = 0, V'(s) = 0, (19) 

where 

4 _ 4 

^(s) = s2(s - /i)2 + (p2 _ k)s^ - 2gs + ^ J . 

Under the condition /3 = {p^ = = a^) we have ip{s) = s(p{s). Thus, relations (fTO)) are similar to 
conditions (fT^ . Therefore, the family of trajectories on the set D* generalizes the family of the especially 
remarkable motions of the 4**^ Appelrot class. 
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3 The equations of integral manifolds 

According to (US]) the finctions S,T form the complete set of the first integrals on D*. In particular, 
the system of equations defining any integral manifold {( E : H{Q = h, K{() ~ k, G{Q = g} is now 
replaced by invariant relations ^ and the equations 

S ^ s, T = T. (20) 

Introduce the complex change of variables [6J generalizing the Kowalevski change for the top in the 
gravity field [3] (i^ = -1) 

xi = (ai - ^2) + i ("2 + Pi), X2 = (ai - /32) - i ("2 + Pi), 

2/1 = ("1 + /?2) + 1(^2 - A), y2 = («! +/32) -i(a2 ^21) 
zi = a3+i^3, 22 = a3-i/33, 
Wi=u;i+iuj2, W2 = tjJi — iuj2, W3 — uj^ . 

In variables ([2T|). system ([6]), ((20)l can be presented in the form 



(?/2 + 2s)wi + a;iW2 + Z1W3 = 0, 

a;2Wi + {yi + 2s)w2 + Z2W3 = 0, , , 

X2Z1W1 + X1Z2W2 + {t - X1X2)W3 = 0, 

2swiW2 — {xiX2 + Z1Z2) + r = 0. 

These equations must be added by geometrical integrals ([2]), which in variables (|2ip can be written as 
follows 

zl + xiy2 = r'^, zl + X2yi^r'^, X1X2 + yiy2 + 2ziZ2 = 2p'^ . (23) 

The space of variables (|2ip has dimension 9 regarding the fact that the following pairs must be 
complex conjugate X2 = STf, 2/2 = yi, Z2 = ^T, and that u;3 is real. Seven relations (|22p . (|23p define then 
the integral manifold. In the case when the the integrals S, T are independent on this manifold it consists 
of two-dimensional tori bearing quasi-periodic motions. 

4 Bifurcation diagram 

Introduce the integral map J of the the dynamical system induced on the closure of the set D*, 

J(C) = (5(C),T(C))eR^ CeCl(D*). 

Due to the obvious compact character of the inverse images of the points of the bifurcation diagram 
E of the map J coincides with the set of its critical values. 

Theorem 3. The bifurcation diagram of the map 

J = 5 X T : C1(D*) ^ (24) 

consists of the following subsets of the (s,T)-plane: 
1°) r = (a -)- 6)2, s e [-a, 0) U [b, +00); 
2°) T = (a - 6)2, s e [-a, -6] U (0, -l-oo); 
3°) s = -a, r ^ (a - b)^; 
4°)s^~b, T^{a-b)^; 
5°)s^b, T^{a + b)^; 
6°) s = a, r s; {a + b)^; 
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'(a2- 


s2)(62 _ s^), 







7°) T = 0, s e (0,+oo); 
8°) r = a2 + fe2 _ 2s2 + 
9°) T = a2 + 52 _ 2s2 - 

10°) T = a2 + 52 _ 2s2 + 2^(s2-a2)(52_;,2)^ g £ +00). 

Proof. Any point of dependence of equations ^ are considered critical for the map (|24p by definition. 
Take the points of trajectories ([5]) that belong to the closure of the set D* and calculate the corresponding 
values of J. These values must then be included in the bifurcation diagram. We obtain those values (s, r) 
for which equations p8| give ([9]), (fTO| . It is easily checked that the half- line pO|) completely belongs to 
the surface given by p^ . It corresponds to the case 1°. Consider the half- line defined by ([9]). The points 
of it lie on (fT8| only if /i2 ^ 4a6. The corresponding set is given by 2°. The segment of ^ in the limits 

-2\foh < h < 2Vab 

is the one-dimensional part of the bifurcation diagram of the map ([5]). For corresponding trajectories ^ 
the value s is not defined. It means that such trajectories are isolated from the set D*. The isolated 
points in the bifurcation diagrams of the reduced systems or of the systems restricted to iso-energetic 
surfaces were met before only in the Clebsch and Lagrange cases. 

To find critical motions in D* use system ((22)) . (|23l) . Introduce the variables x, z, 



x^ = X1X2, = z\Z2. (25) 
It follows from the last equation that 

2/1Z/2 = 2/ -x" - 2z^ (26) 

and the first two equations give 



(zi + Z2f 2r2 - (xiy2 + X2yi) + 2z2, 
{zi - Z2f = 2r2 - [xxy2 + X2yi) - 2z'^. 



(27) 



Eliminate zi , Z2 in 
Then using we obtain 

Denote 

From (1271), dMl) we have 



(r^ - Xiy2)ir^ - X2iji) = z'^ . 
r^{xiy2 + X2yi) ^r'^ + 2p^x^ - {x^ + z^ f. (28) 



^±{x,z) = (a;2 + z2±r2)2-2(p2±r2)a;2 



r2(zi +Z2f = '^+{x,z), r^{zi - Z2f = <^-{x,z). 
Therefore, x, z satisfy the inequalities 

$+(a;,z)^0, <^-{x,z)^{). (29) 

Notice that the equilibria of system ^ are included in the family of motions Vl. In all other cases 
the determinant of the first three equations \22\ in Wi {i = 1, 2, 3) is identically zero. Eliminating , z|, 
and 2/1 j/2 with the help of ([23| . ([26]) . we obtain 

2s[{r'^xi - ryi) + {r^X2 - rya)] = -r'^{xiy2 + X2yi)+ .„„s 
+2[2s2(r - a;2) + p^{t + x^) ~ t{x^ + zJ)]. ^'^^> 
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On the other hand, and (1261) yield 

(r^xi - Tyi){r^X2 - Ty2) = r^x^ + r(2p^ - - 2z^) - r'^T{xiy2 + a;22/i). (31) 

Denote 

Then from the second relation (fTSl) the identity follows 

Introduce the complex conjugate pair 

^^^r^xi-ryi, = r'^X2 - Ty2. 

Eliminating the expression xiy2 + X2yi in (|30ll . (|3T|) with the help of (|28ll . we obtain the system 

2s(mi + M2) = (a;^ - rf - 4s^x^ -a + As'^t, , , 

MiA^2 = T(a;^ + — t)^ + ctx^ — tct. ^ ' 

Choose 



Ai = \/2s^i + a, A2 = \/2s[i2 + cr 
to be complex conjugate. Then system (|33l) takes the form 

(Ai + \2f = z), (Ai - Az)' = z), 

where 

*±(x, z) = (x^ + - r ± 2sx)^ - 4s^x^ 

It is solvable if 

^+ix,z)^0, ^^ix,z)^0. (34) 

The system of inequalities (P^ . (IMl) defines the region of possible motion (the RPM) in the (a;, z)- 
plane. The RPM is the projection of the integral manifold Jg^r- For given s,t the initial phase variables 
are algebraically expressed in terms of x, z. The bifurcation diagram corresponds to the cases when the 
RPM undertakes qualitative transformations as its parameters s,r change. 

Introduce the local coordinates si, S2 in {x, z)-plane: 

x"^ + z^ + 



Sl = , S2 



2x 2x 
Inequalities (^H) are immediately solved 

sl ^ a^, sl ^ b^. (35) 

The corresponding region in the (x, z)-plane is shown in Fig. 1 for the first quadrant. We also point out 
the coordinate net (si,S2). 

Let Hi be the rectangle in (si, S2)-plane with the vertices si = ±a, S2 — ±b. To solve system (IMl) 
express 



2^ - T = [Si + S2 j(si - S2)]x, 



where 



*+(a;, z) = a;^A+A_, *_(a;, z) = x'^M+M-, 
T — 2sx 

A±(si, S2) = Sl + S2 2 — ~ ^2) ± 2s, 

T + 2sx 

M±(si, S2) = Sl + S2 5 — (si - S2) ± 2s. 



T 

I Q-) — ^ I'r 

(36) 
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a-b 



a + b 



Figure 1: The admissible points in the [x, z)-plane. 



It foUows from JM]), ^ that 



A+(si,S2)A_(si,S2) ^ 0, M+(si,S2)M_(si,S2) 0. 



(37) 



Consider the paraUelogram 112 bounded by the Hues A± = 0, M± = 0. The solutions of system ([57]) 
fill two half-strip regions starting at the sides of 112 belonging to the lines A± = 0. The example of the 
RPM in the (si, S2)-plane, i.e., the set of solutions of inequalities (|35p. ([57)) is shown in Fig. 2. 

The further investigation is purely technical. The bifurcations of the RPM's take place in one of 
the following cases: the vertex of one parallelogram out of 111,112 resides on the boundary of another; 
the sides of the parallelograms Hi , 112 happen to be respectively parallel (the vertices of the RPM go 
to the infinity); the half stripe region degenerates and becomes a half line. Finding all such cases gives 
the equations for s, r pointed out in the theorem. Let A denote the set defined by these equations in 
R^(s, r). Considering the connected components of R^(s, t)\A we ignore those of them which correspond 
to the empty RPM's. The rest of components (the admissible regions in the integral constants space) are 
shadowed in Fig. 3. The bifurcation diagram consists of those segments of A that are boundaries of the 
admissible regions excluding the parts of the axis s — Q since this value of s is not admissible by virtue 
of . This way we obtain the inequalities needed. The theorem is proved. □ 



Due to the second equation in (j33p each trajectory is represented by some curve on this surface. 

Obviously, the constants t, a cannot be simultaneously negative. Therefore, surface (|38p has two 
families of rectilinear generators. 

The introduced constants satisfy the following two identities 



5 On the possibility of the separation of variables 



Denote 




2 /-2 , 2 



(38) 



(39) 
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Figure 2: The example of a region of possible motions (a = 1, b = 0.4, r = 1.2, s = —0.6). 



It is easily seen that by virtue of ([5^ . the equations <i>± = 0, = in the (a:,^)-plane define 
the family of lines tangent to the cross section of surface ([55)1 by the plane ^ — 0. Such line is then the 
projection of some generator of surface (|38p . Since each point on surface ([55]) belongs exactly to two 
generators the parameters of the latter can be chosen as local coordinates in the region in the {x, ^)-plane 
covered by surface p8|) . 

Not regarding any reality conditions, put formally 

— UV + 1 /—U — V 

i^^a ■ , x = Vt — ■ — . 

u-\- V U + V 

After some simple transformations we get 

^+ = 1 — ^-T^'Pi{u)'Pi{v), = 7 — ^-T^V2(w)v52(w), 

[u + vy [u + v)^ 

(u + vy [u + vy 

ipi{w) = + + 2(t + r2)w, ip2{w) = ^/^{l + w^) + 2(t - r^)^, 

■01 (w) — ^/(t{\ + w^) + Asxw, ip2iw) = + w'^) — Asxw. 

In the plane of the variables u,v inequalities ([^5]) . ([M]) define the set of rectangles with the sides 
parallel to the coordinate axes. The fact that each connected component of any integral manifold is 
represented by such a rectangle (and in the case of a bifurcation by a segment or a pair of rectangles 
having a common side) means that in these variables the equations of motion must separate. The 
corresponding calculations are too long for the restricted volume of this article and will be presented in 
another publication. We only point out the connections with the above results. 



where 
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Figure 3: Region of existence of motions in the (s,T)-plane. 



Consider the polynomial 

Q(w) = (y3i(w)(p2(w)V'l(w)V'2(w) (40) 

and find all the cases when it has a multiple root. The resultant of Q(w) and Q'[w) in w is (up to the 
constant multiplier) 

sVi2(r2 - 2/r + r4)i4[2s2 - (j? - r2)]4[2s2 - (j? + r'')\\T^ - 2(j? - Is'^Y + r^f. (41) 

As it was already mentioned, by virtue of equations (|18p at the considered family of motions we 
have s ^ 0. The rest of cases when expression (|4T|) vanishes lead to the equations listed in Theorem 3. 
Therefore, the bifurcation diagram found above is the part of the discriminant set of polynomial (|40p . 
Such phenomenon is also typical for the systems with algebraically separating variables. 
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